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ABSTRACT 



By adapting the streamline model developed by Caudle and LeBlanc 
for flow through porous media, a streamline model for in-situ combustion 
was developed. 

Initial tests were made to compare the streamline model for in-situ 
combustion with actual field trials conducted by oil companies. The results 
were extremely close. As a result, in-situ combustion experiments may 
be examined and optimized to obtain the best possible results prior to im- 
plementation in the field. 

Conducting numerous tests by changing assumed values indicated 
that areas in the process must be re-examined in order to determine the 
economic feasibility of using in-situ combustion to recover the residual 
oil remaining in the nation's reservoirs, kerogen from oil shale, and the 
hydrocarbons from tar sands. 



iv 



TABLE OF CONTENTS 



PAGE 

ACKNOWLEDGMENTS iii 

ABSTRACT iv 

TABLE OF CONTENTS v 

LIST OF TABLES vi 

LIST OF FIGURES vii 

CHAPTER 

I INTRODUCTION 1 

n THE STREAMLINE MODEL 3 

m IN-SITU COMBUSTION MODEL 9 

IV EXPERIMENTAL INVESTIGATION 26 

V CONCLUSIONS AND RECOMMENDATIONS 54 

APPENDICES 

A FLOWCHART 57 

B COMPUTER PROGRAM 62 

C NOMENCLATURE 68 

BIBLIOGRAPHY 74 

VITA 



V 



LIST OF TABLES 



TABLE PAGE 

I Recovery Data for Basic Streamline Model of 

In-Situ Combustion in an Unbounded System 30 

n Recovery Data for Streamline Model Utilizing 

Moss Relationship Between Advancement of 
Combustion Front and Injected Air Flux 33 

m Recovery Data for Basic Streamline Model of 

In-Situ Combustion in a Bounded System 46 

IV Characteristics of Sun Oil Company Field Test 48 

V Characteristics of Magnolia Petroleum Company 

Field Test 50 



vi 



LIST OF FIGURES 



FIGU RE PAGE 

1 Unbound Five Spot Pattern Showing Generated 

Streamlines from Wells Having Equal Rates 7 

2 Typical Streamlines in a Two Well System 8 

3 Typical Flow Tube 8 

4 Gate’s View of the In-Situ Combustion Process 10 

5 Allred’s View of the In-Situ Combustion Process 11 

6 In-Situ Combustion as Used in the Modified 

Streamline Model 12 

7 Effect of Injected Air Flux on the Combustion 

Front Rate of Advance 14 

8 Process Zones in a Non-Linear System 21 

9 Process Zones in a Linear System. 21 

10 Streamlines in the Basic Reservoir Used for In- 

Situ Combustion 28 

11 Location of Oil Zone in Basic Test at Breakthrough 31 

12 Combustion Front Location Using Martin and Moss 

Relationship Between Injected Air Flux and Com- 
bustion Front Rate of Advance 32 

13 Effect of Variation in Hydrocarbon Consumption 35 

14 Effect of Variation of Zone One Mobility 37 

15 Effect on Oil Production by Variation in Injected 

Air Flux 39 

16 Effect on Time of Breakthrough by Variation in 

Injected Air Flux 40 

17 Combustion Front Location at Breakthrough with 

Variation in Injected Air Flux 42 



LIST OF FIGURES (Continued) 



FIGURE PAGE 

18 Typical Pattern Showing Imaging Technique 43 

19 Location of Combustion Front in a Bounded 

System 45 

20 Streamline Model Burned Zone versus Sun Oil 

Company Field Test 49 

21 Streamline Model Burned Zone versus Magnolia 

Petroleum Company Field Test 51 



viii 



CHAPTER I 



INTRODUCTION 



In 1972 the United States became aware of an energy crisis. The 
crisis took the form of shortages of fuel for production of electricity and 
heating fuels for homes and businesses. The causes are many and varied, 
but may be traced to one reason — economics. Since the start of oil pro- 
duction, the cost per barrel has been kept so low by artificial barriers that 
the return of the dollar invested in oil or gas exploration and production is 
below the return of investments in other areas. Thus the energy crisis is 
not the lack of fuel per se, but the lack of economic fuels. 

All persons intimately related with the oil industry know that existing 
means of production may optimistically recover fifty percent of the oil in 
the nation’s reservoirs. Thus one could conservatively estimate that over 
90 billion barrels of oil still remain in established reservoirs, since the 
total production of the United States from 1920 to 1969 was just over the 90 
billion barrel mark (1). In addition, potentially 1,781 billion barrels of 
kerogen from oil shale and 399 billion barrels of oil from tar sands exist in 
the United States (2). Based on an estimated rate of consumption of 46 mil- 
lion barrels per day by 1980, the above resources could supply the United 
States for over 135 years. 

The fallacy with the above conclusion is that the price of a barrel 
of oil would have to rise drastically to make production of the residual oil, 
oil shale, and tar sands profitable if present production concepts are main- 
tained. For example, present methods of oil production from oil shale is 
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by mining the oil shale and retorting the shale above ground. Not only is 
the process environmentally distasteful, but the resulting cost per barrel 
of oil is between four and five dollars (2). 

During the past ten to fifteen years, one process has received exten- 
sive investigation — in-situ combustion, often called fire flooding or under- 
ground retorting. Through in-situ combustion, residual oil in existing oil 
fields could be recovered. Synthetic oil from oil shale and tar sands could 
be recovered without strip mining, and possibly at lower costs than the 
present. 

To the author's knowledge no attempt has been made to simulate in- 
situ combustion by use of the Caudle -LeBlanc Streamline Model (3). Ex- 
pressed simply, the Caudle-LeBlanc model is a mathematical representation 
of the movement of one typical particle of oil. The following study describes 
the basis of the model and the modification of the model to include the specific 
characteristics of in-situ combustion, with the concluding sections comparing 
the model results with the results of actual field tests conducted by various 
organizations. 



CHAPTER II 



THE STREAMLINE MODEL 



The theory and development of the streamline model has been re- 
viewed by many different authors (3, 4). For this reason only a brief review 
of the model will be shown. Following this brief introduction the model 
modifications will be shown. 

The Basic Streamline Model 

The basic equations were developed from Darcy’s Law and included 
the following assumptions : 

1. negligible gravitational effects inside the reservoir 

2. incompressible homogeneous fluids 

3. horizontal system, infinite in extent 

4. homogeneous, isotropic media 

The most general form of Darcy's Law is as follows: 



kA dp 
~ ~d2 



( 2 . 1 ) 



It is also known that the fluid flux is equal to the flow rate divided by the 
area through which it flows, thus Equation 2. 1 may be written as 



S_ - - iL 

~A ~ IJL dC 



( 2 . 2 ) 
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In order to determine the streamline path it is necessary to know its velocity 
in the x, y, and z directions. Since it is known that the velocity is the fluid 
flux divided by the porosity of the material through which the fluid is flowing, 
Equation 2. 2 may be written as: 



Velocity = v = — ^ (2.3) 

0 0 /i d^ 

The second important equation is the continuity equation. It is a partial 
differential equation derived from the material balance around a point in 
three-dimensional space. This equation is written as follows: 



a (pux) a (puy) a (pu^) 

+ + = ~ 0 "jr" 

ax ay az dt 



(2.4) 



When the volumetric flux components of Darcy's Law, Equation 2.2 in the 
X, y, z directions, are substituted in Equation 2.4 the following equation 
results: 



3 

ax I p ax 



^ a 

9y \ M 9y 




' pk^ 3^ ^ 
^ p az j 




(2.5) 



Applying the assumptions made earlier. Equation 2. 5 becomes: 



a^p ^ ^ 

ay2 



0 



( 2 . 6 ) 



which has become known as LaPlace's Equation. By transforming to cylin- 
drical coordinates and allowing radial flow around a point, the equation may 
be solved for the flow potential in a horizontal, isotropic, homogeneous 
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porous medium. Restoring the resulting equation to rectangular coordinates 
and allowing for a multiwell system, Equation 2.6 becomes: 

n 

p(x, y) = P„ - Qi tn [(x-Xj)^ + (y-y/J ( 2 . 7) 



where is the mean reservoir pressure. 

By differentiating Equation 2.7 with respect to x and substituting the 
resulting derivative into Equation 2. 3 the velocity at any point in a rectan- 
gular system becomes 



(x, y) = 



n 

I 



q (x-xj) 



l = l(x-Xj)2 + (y-y.) = 



( 2 . 8 ) 



In a similar manner, the velocity in the y direction becomes 



n 



Vy (x, y) ^ — y 

y o.TT thYi . 



q(y-yi) 



i = l(x-Xj )2 + (y-yi)2 



(2.9) 



In order to follow the particle along a particular streamline, the dis- 
tance of travel must be known. Given a finite amount of time, the distance 
a particle travels at a constant velocity may readily be calculated as 



ds = vdt or As = vAt (2.10) 

Assuming a particle was initially located at point (x^ , y^) , the location of 
the same particle after time At could be calculated as 
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’‘i+1 = 

^i+i = yi * 



( 2 . 11 ) 



By starting on the wellbore of the injection well and systematically 
reiterating Equations 2. 11, a trace of any particle's movement may be made. 
Figure 1 shows the paths of particles, from a single five -spot pattern with a 
central injection well. The time to reach a specific point may be determined 
by the summation of the At's required to cause sufficient movement of the 
particle to reach that specific point. 

In a given two-dimensional fluid system, an infinite number of par- 
ticle paths (streamlines) could be used to describe the flow from a source to 
a sink. However, to simplify computation, a representative number of 
streamlines may be chosen to represent all fluids emanating from the source. 
A few of the streamlines from a source are shown on Figure 2 (5). By 
dividing the distance between two adjacent streamlines into two equal parts, 
the streamline may be said to represent the flow of all particles within the 
dividing lines. A streamtube has thus been defined and is shown on Figure 3. 
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O Production Well 
• Injection Well 



FIGURE 1 Unbounded Five-Spot Pattern Shwoing Generated 
Streamlines from Wells Having Equal Rates 
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FIGURE 2 Typical Streamlines in a Two-Well System 




FIGURE 3 Typical Flow Tube 



CHAPTER m 



IN -SITU COMBUSTION MODEL 



Prior to developing the modifications necessary to convert the stream- 
line model to a streamline model for in-situ combustion, a discussion of the 
process of in-situ combustion is warranted. 

In-situ combustion is not a new process. Recovery by a combustion 
or heat wave process was patented in 1923. F. A. Howard was granted a 
patent on a process in which air and a combustible gas were pumped into an 
injection well and ignited (6). The method involves ignition of the formation 
in an injection well, followed by propagation of a combustion front through 
the reservoir. Combustion is maintained by the injection of an oxygen- 
containing gas, such as air, to react with reservoir hydrocarbons. As the 
flame front, which is really glowing embers of the reservoir material, pro- 
gresses through the reservoir, oil and formation water are vaporized, driven 
forward in the gaseous phase and recondensed in the cooler part of the for- 
mation. In turn the condensed fluids push oil into the producing wells (7). 

Since 1923 field experiments have been conducted by the Russians in 
1935 (8) and others in Oklahoma (9) and Kansas (7), with the most recent 
test occurring in the Green River Oil Shale deposits by the U.S. Bureau of 
Mines (11, 12). The objective of the early tests and those occurring in the 
later years was to recover the residual oil left after primary drive and 
secondary methods had been utilized. In the 1950 ’s in-situ combustion ex- 
periments were conducted on oil shales with the intent of retorting the 
shale underground as an alternative to mining the shale and retorting above 
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FIGURE 4 Gate’s View of the In-Situ Combustion Process^ 



INJECTION PRODUCTION 
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FIGURE 5 Allred’s View of the In-Situ Combustion Process 
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FIGURE 6 In-Situ Combustion as Used in the Modified 
Streamline Model 
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ground. Other tests were conducted in highly viscous oils with the intention 
of reducing the viscosity of the oil in order to make the oil more mobile and 
thus producible in a smaller period of time (10). 

The complex nature of the in- situ combustion process, with it’s 
many inter-related displacement mechanisms, makes exact interpretation 
of field data impossible. This is understandable considering that even two- 
phase flow of fluids through a porous homogeneous medium, which has been 
studied for many years, cannot be adequately described. As a result, as- 
sumptions must be made based on e?q)erience and the results of many trials. 
Gates and Ramey (10) conclude that, as a result of their interpretation of 
field tests in California, Figure 4 represents the in-situ combustion process. 
Allred’s (13) beliefs are somewhat different and describe the in-situ com- 
bustion process as shown on Figure 5. With the exception of the gases 
shown in Figure 5, both representations are the same. This author has 
simplified the two figures to enable utilization in the streamline model simu- 
lation. The in-situ combustion process as utilized in the simulation is 
shown on Figure 6. 

Along with the composition of the various zones, the velocity of the 
fronts is important. Since the burn front is the driving mechanism in the 
process, its velocity is the basis for determining the velocities of the other 
zones. Two experiments indicate that the injected air or gas must reach a 
certain rate before movement of the front will begin. Figure 7 represents 
the results of the aforementioned ejqseriments (9, 14). Further experiments 
may show that the rate of advancement of the combustion front depends on 
the oil content of the reservoir, reservoir porosity or some other variables. 
What is important at this time is that the relation between air flux across 
the burn and the rate of movement of the burn does exist. 
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FIGURE 7 Effect of Injected Air Flux on the Combustion Front Rate 
of Advance 
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Other such relations and assumptions had to be established in the 
modification of the streamline model, but to ensure their complete compre- 
hension they will be discussed as they affect the streamline model modifi- 
cation. 

In developing the original streamline model we assumed that the fluid 
was incompressible. The effect of this assumption is extremely close to reality. 
Merrick (18) compared the results of flow of real gases, ideal gases, and 
incompressible fluids through a porous media and determined that the dif- 
ferences in results were extremely small. This assumption will also be 
made for the in- situ combustion model of flow in a reseiwoir, but to utilize 
the relationship of the advancement of the combustion front and the injected 
air flux, the injected air flux must be in surface volumes. To properly 
change the injection rate to surface volumes we must look at the gas laws. 

For an ideal gas 



where V is the volume of a given number of gas moles n at absolute pres- 
sure p and temperature T. R is the universal gas constant. For a non- 
ideal gas the relation is 



where z is the compressibility factor of the gas. The compressibility is 
a function of the pressure in the reservoir Pj^. Thus, Equation 2,3 de- 
scribing an ideal gas becomes 



pV = nRT 



(3.1) 



pV = znRT 



(3.2) 




V 



(3.3) 
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where is the compressibility factor of the fluid injected into the reser- 
voir, and is the average reservoir pressure in atmospheres (standard 
pressure is one atmosphere) and the surface and reservoir temperature are 
equal. 

Equations 2.7, 2.8, and 2.9 must also be modified since they 
were developed for noncompressible flux. With an ideal fluid, the volume 
at standard conditions equals the volume at reservoir conditions. Since we 
are using a compressible fluid, the quantity at surface conditions becomes 



q 



sc 




^res 



(3.4) 



Since it is also known that the mobility of a fluid is equal to its permeability 
divided by its viscosity. Equations 2.7, 2.8, and 2.9 become 



P 



'm 



^m 



4nhp„jXs 





(3.5) 



n 



Vs, 






^ (x-xi) 

I^scj - ^2 



27Th(^Pm i=i (x-x^) + (y-yj)' 



(3.6) 



n 



Vs 



"m 






(y-yi) 



y 2;rh0Pjjj i (x-Xj) + (y-yj)' 



(3.7) 



From the above equations the pressure and velocity of the standard fluid may 
be determined. The standard fluid in this thesis will be the air injected into 
the reservoir. It is called the standard fluid to keep the process as general 
as possible . In addition, the standard fluid has particular significance in the 
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conductivity ratio which will be explained in the coming pages. 

To determine the amount of production, the velocity of the burn is 
necessary since, as will be shown later, the movement of all the zones 
containing oil are dependent on the movement of the burn. Equations 3.6 
and 3.7 will generate the velocity of the standard fluid, which in this case 
will be the air injected into the reservoir. Knowing the velocity of the air, 
and by the relationship of the air velocity and burn front velocity as expres- 
sed in Figure 7, the velocity of the burn may be computed. From Figure 7 
the velocity of the burn is 





= SL Vg + C 


(3.8) 


x 


X 






= SL Vg + C 


(3.9) 



y y 

where SL is the slope of the line and C is the line constant. Note that 
at approximately thirty feet per hour, the relation between the burn and the 
standard fluid changes. At this point the values of SL and C will be 
changed. In the streamline computer program these values are used as input 
data to provide as much flexibility as possible. As previously mentioned, 
insufficient testing has been completed to prepare exact figures for the burn 
front movement with regard to the rate of injection of air, but the fact that 
the relation does exist must be considered. 

The initial streamline model developed in Chapter II was for a single 
fluid. For a multifluid condition, such as a waterflood or in-situ combustion 
where different fluids exist, the concept of the conductivity ratio must be 
applied. Caudle (15) presents a detailed explanation of the concept and thus 
it will not be discussed here. The conductivity ratio is determined as fol- 
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lows: 



y 




(3. 10) 



where y is the conductivity ratio, is the pressure at the source or in- 
jection well, Pp is the pressure at the sink or production well, Xg is the 
mobility of the standard fluid, and Xj^ is the mobility of the fluid in the k 
zone which has a pressure pj^ at its interface with the zone preceding the 
k zone. L is the number of different zones. 

Substituting Equation 3. 3 in Equations 3. 8 and 3. 9 and multiplying by 
the conductivity ratio 



As a result the velocities (x and y directions) may be computed for 
the burn front. The velocities of the other three fronts will be similar, but 
will differ by some factor due to the accumulation of oil in front of the burn. 
By referring to Figure 8 it may be seen that the velocity across point B 
will yield a quantity of a known amount if the cross-sectional area of 
point B is known. By the same rationale a fluid moving across point 2 
will yield a quantity of fluid Q 2 of known amount if the area and velocity 
of point 2 are known. Additionally, if a constant rate of injection is estab- 




(3.11) 



y 




(3. 12) 



lished, the quantity of fluid passing point B will, in some way, be related 
to the quantity of fluid passing point 2, Mathematically, 



% = ''b^ 


(3. 13) 


Q2 = ''2^2 


(3.14) 


*^2 = ^2% 


(3. 15) 


Thus, by substituting Equation 3. 13 in Equation 3. 15 




«2 = ^2%^ 


(3. 16) 


Considering the constant rate of injection, the following conclusions may be 
made : 


= ''Sb\ 


(3.17) 


% = ''S2*2 


(3. 18) 



Equating Equations 3. 17 and 3. 18 and solving for A 2 , 

''Sb^b 

A 2 = (3.19) 



Substituting Equation 3. 19 in Equation 3. 16, 
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Solving for the velocity at point 2, 




(3.21) 



By similar reasoning point 3 on Figure 8 would have a velocity of 




(3.22) 



As shown on Figure 6, the fourth zone is composed entirely of gas, and by 
a piston-like displacement the zone is pushed by the air injected into the 
reservoir. Thus the velocity of zone 4 is that of the standard fluid at zone 
4 as computed by Equations 3.6 and 3.7. 

Considering all the x and y components of the velocities of zones 
2, 3 the movement of these two points may be computed as follows: 



X 




(3.23) 



y 




(3.24) 



X 




(3.25) 
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FIGURE 8 Process Zones in a Non-Linear System 




FIGURE 9 Process Zones in a Linear System 
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V 




^ 3 % 





(3.26) 



To compute the values of F 2 and F^, a movement of the burn must 
be considered. Referring to Figure 9, a movement of the burn from B to 
B’ will cause the displacement of all oil existing in the area between B and 
B’ that is not consumed as fuel in the combustion process. Thus 



Oil moved 






(1. 0 - Sg^ - 

(1.0 -Sg.) 



(3.27) 



Since 



q = 






dt 



(1.0 -Sg.) 



(3.28) 






^i </>eff 



(3.29) 



dVg I " O (1-0 - fij) (1. 0 - Sg^ - S^y^ 



dt A 0 



eff 



(1-Sgp 



(3.30) 



Since 



dVg/Aj 

dt 



(3.31) 
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Equation 3,30 may be written as 



■^eff (1-0 -Sgj) 

where 0^^^ is the effective porosity through which the fluid must move and 
is computed as 



Vj,. (3.32) 

1 



<’et£ = '^g(l-«-Sg-S„) (3.33) 

As was discussed earlier, the velocity of the burn is proportional to the 
velocity of the air flux across the burn front. Thus if we consider point 2 
as the i-th point, 

(3.34) 

(3.35) 

and Equation 3.32 becomes 




(0T ~ ( 1 • 0 - f|j) ( 1 - Sg^ ~ S^^) 

0eff(l-Sg2) 



V V, 

S2 ''b 



(3.36) 



Comparing Equation 3. 36 with 3.21, we find that F 2 is 
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(«^^-0g)(l.O-fb)(l.O-Sg^-S^2) 



(3.37) 



By mutual displacement, the movement of point 2 to 2' has moved 
point 3 to 3' by an amount equal to F 2 . Thus the movement of point 3 to 
3’ is 



V 



3' 



^s3’ 

^sb 



(3.38) 



The movement of point 3’ to 3" caused by the movement of the burn and 
not the movement of point 2 to 2' is given by Equation 3. 32 as 



^3' 



Combining Equation 
to 3” , 






S 3 , b 



•^eff^^*® ‘^go,^ 



(3.39) 

3. 38 and 3. 39 to obtain the total movement of point 3 



(0,p - 0g) ( 1 . 0 fj^) ( 1. 0 - Sgg - 3^2^ 
0 ff(l-S„ ) 

^eff^ gg' 



+ F, 



Vc 

S 3 b 



V 



% 



(3.40) 



Comparing Equation 3.40 to Equation 3. 22, F^ is found to be 
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+ F. 



(3.41) 



Substituting the equations developed in this chapter into the streamline 
model developed by Caudle and LeBlanc (3) will result in an in-situ combus- 
tion streamline model. 



CHAPTER IV 



EXPERIMENTAL INVESTIGATION 



Sensitivity to Assumptions 

As was discussed earlier, little is known about the in- situ combustion 
process. Present technology prevents us from obtaining cores from the 
reservoir as the process moves through it. As a result, numerous assump- 
tions had to be made. The first assumption is the relationship between the 
advancement of the combustion front and the injected air flux. The effect on 
production and time of breakthrough by changes in this relationship must be 
determined. The second assumption is the amount of reservoir hydrocarbon 
consumed in the propagation of the combustion front. As may be seen in 
Equations 3.36 and 4.40, the velocity of zones 2 and 3 are directly re- 
lated to the amount of reservoir hydrocarbon consumed. The third variable 
is the mobility of the various zones within the process. 

In order to determine the sensitivity of the streamline model to vari- 
ations in these assumptions, a basic test was conducted with the following 
data. 



1. The rate of advancement of the combustion front in relation to 
air flux is in accordance with the curve in Figure 7 established by Martin 
(14) because it shows a realistic slowing of the combustion front as the air 
flux declines. 

2. The saturation of gas, oil, and water in each zone is 
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Zone 1: Sg =1.00 

Zone 2: Sg =0,33 

So = 0.34 

s„ = 0.33 

Zone 3: Sg =0.50 

So = 0.30 

S^ = 0.20 

Zone 4: Sg = 1.00 



3. The rate of injection was established at 600scf/hr. 

4. The mobility of each zone was assumed to be 25, 11, 12, 25, 25 
darcy/centipoise for Zones 1, 2, 3, 4, and5 respectively. The wisdom of 
using these mobility ratios in light of the saturation and make-up of each 
zone may be dubious; however, tests were conducted to determine the effect 
of changes in zone one mobility only. The others were not examined due to 
reasons to be explained later. 

5. The amount of oil consumed in the propagation of the combustion 
front was assumed to be ten percent of the total hydrocarbons in the reser- 
voir. 

6. The average amount of oil in the reservoir was assumed to be 
27 gallons of oil per ton of reservoir material, 

7. Pattern utilized was an unbounded inverted five spot, 100 feet 

square . 



The above data were incorporated into the computer program shown 
in Appendix B. The plot of the streamlines is shown in Figure 10, and the 
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FIGURE 10 Streamlines in Basic Reservoir 
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production data are shown in Table I with the location of the oil bank at 
time of breakthrough in Figure 11. In subsequent tests the location of the 
combustion front will be compared with the basic model plot in order to 
show the amount of the reservoir swept by the burn. 

The relationship between the advancement of the combustion front 
and the air flux is fundamental to successful simulation of in-situ combus- 
tion. The two curves showing this relationship in Figure 7 have been pub- 
lished as a result of experimental studies (9, 14). As a result, these curves 
must be assumed to include such effects as oxidation potential of the oil, 
heat conduction and heat convection. The basic program was computed 
using the relation developed by Martin (14). To determine the significance 
of the difference in the Martin and Moss (9) curves, the second test was to 
use all the data from the basic test with the Moss curve. Figure 12 com- 
pares the results of the tests, where the basic data using the Martin curve 
is a solid line and the second test using the Moss curve is a broken line. 
From the start of ignition to 129 hours, the fronts are propagated radially. 
The positions of the fronts shown in the upper right quadrant is the position 
of the combustion fronts at breakthrough. The number of hours shown in 
that quadrant is the time of breakthrough of the indicated front. Table II 
shows the production data for the second test. The difference in the break- 
though time and the total production may be explained by review of Figure 7. 
At high air flux rates, the front of the basic test moves at a faster rate. 
Production data confirms this in that the initial oil production is 97 percent 
greater in the basic test, indicating a larger area swept in the same amount 
of time. As the air flux rate drops due to the expanding surface area of the 
burn front, the Martin curve declines more rapidly than the Moss curve. 
After 629 hours the combustion front of the basic test has reached the pro- 
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TABLE I 



RECOVERY PERFORMANCE FOR BASIC TEST 



Gas Produced 
(ft^) 


Oil Produced 
(bbls) 


Water Produced 
(bbls) 


Time 

(hrs) 


0.0 


0.0 


0.0 


4.2 


0.0 


0.0 


0.0 


45.2 


0.0 


0.0 


0.0 


87.5 


53.3 


3.3 


2.4 


129.0 


143.0 


10.0 


8.3 


171.0 


252.0 


18.2 


15.8 


212.0 


376.0 


27.8 


24.7 


254.0 


506.0 


38.2 


34.4 


296.0 


649.0 


49.4 


44.8 


337.0 


779.0 


59.6 


54.6 


379.0 


948.0 


72.6 


66.4 


421.0 


1070.0 


83.0 


76.4 


462.0 


1190.0 


92.9 


86.0 


504.0 


1360.0 


106.0 


98.0 


546.0 


1480.0 


115.0 


106.0 


587.0 


1580.0 


122.0 


114.0 


629.0 
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FIGURE 11 Location of Oil Zone in Basic Test at Break- 
through 



Basic Test (Martin Relationship) 
Second Test (Moss Relationship) 




FIGURE 12 Combustion Front Location Using Martin 
and Moss Relationship Between Injected 
Air Flux and Rate of Frontal Advance (14, 9) 



TABLE II 



RECOVERY DATA FOR SECOND TEST 



Gas Produced 

(ft5) 


Oil Produced 
(bbls) 


Water Produced 
(bbls) 


Time 

(hrs) 


0.0 


0.0 


0.0 


4.2 


0.0 


0.0 


0.0 


45.8 


0.0 


0.0 


0.0 


87.5 


56.8 


0. 1 


0. 1 


129.0 


145.0 


2.8 


1.9 


171.0 


242.0 


7.8 


5.2 


212.0 


356.0 


14.2 


9.4 


254.0 


481.0 


21.8 


14.7 


296.0 


606.0 


30.5 


21.6 


337.0 


734.0 


40.4 


30. 1 


379.0 


864.0 


50 . 6 


39.0 


421.0 


1000.0 


61.7 


49.0 


462.0 


1140.0 


73. 1 


59.4 


504.0 


1280.0 


85.3 


70.6 


546.0 


1430.0 


98.2 


82. 5 


587.0 


1580.0 


112.0 


95.3 


629.0 


1740.0 


126.0 


108.0 


671.0 


19 10 . 0 


140.0 


121.0 


712.0 


2070.0 


154.0 


134.0 


754.0 


2240.0 


169.0 


149.0 


796.0 


2400.0 


184.0 


163.0 


837.0 


2550.0 


198.0 


177.0 


879.0 


2670.0 


211.0 


189.0 


921.0 


2800.0 


224.0 


202.0 


962.0 


2920.0 


235.0 


213.0 


1000.0 
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duction well and the streamlines midway between two production wells have 
stopped. The reason for these streamlines stopping may be found in Equa- 
tion 3. 10 which has been expanded below to include all zones considered in 
this thesis. 



y = — 

X r^ii^ 

S Xg X2 X3 X4 X5 

(4.1) 

As the streamlines between two production wells move out from the injection 
well, the pressure at those points decreases which in turn decreases the 
conductivity ratio. Since the conductivity ratio is directly proportional to 
the velocity of the burn, as shown in Equation 3.11, a decrease in the con- 
ductivity ratio will cause the combustion front to slow and eventually stop. 
The same would be true for the streamlines reaching the production wells, 
except that, as each zone is produced, its term in the denominator of Equa- 
tion 4.1 is dropped thus causing an increase in the conductivity ratio. The 
slowing and stopping of the streamlines between the production wells causes 
the oil front ahead of the burn to stop before reaching the production wells. 

The next area of investigation was to determine the effect of varia- 
tions in the amount of hydrocarbon consumed in the propagation of the com- 
bustion front. Equations 3. 38 and 3.42 show the relationship between the 
velocity of the different fronts to the amount of reservoir hydrocarbon con- 
sumed. As the amount of reservoir hydrocarbons consumed is increased, 
the velocity of each front will be decreased. Figure 13 shows that as the 
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Total Production, bbls 
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FIGURE 13 Effect of Variation in Hydrocarbon Consumption 



Time for Breakthrough of Combustion Front, days 
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amount of hydrocarbon consumed is increased from 10 percent to 60 percent 
the amount of production and the time required for breakthrough increase. 
When the consumption rate is set at 70 percent of the reservoir hydrocarbon, 
the total production and time of breakthrough decrease. 

The results indicate that, as the velocity of the combustion front de- 
creases, the amount of production also decreases. Confirmation of this 
effect should be obtained when the injection rates are varied. However, 
prior to studying the effect of variation in injection rates, the last assump- 
tion should be discussed. 

The third assumption made in the utilization of the streamline model 
is the mobility of zone 1. As the combustion front moves across the reser- 
voir, the water and oil have been driven forward or been used as fuel for 
the combustion process. In shale the clay would shrink much as a brick does 
during firing. Logically, the permeability of the zone behind the combustion 
front should increase and thereby enable the injected gas to pass through the 
zone more easily. The mobility of this zone is 

X = -!L (4.2) 

where k is the permeability of the material through which a gas or liquid of 
a specific viscosity (/r) is passing. With an increase in permeability, the 
mobility should also increase, assuming the viscosity of the gas or liquid 
does not change. The mobility of the first zone has the greatest effect on the 
velocity of the combustion front since it is the first zone which contains the 
standard fluid. Equation 3. 10 shows the rationale for the above statement. 
Figure 14 shows the results for four trials using a mobility of the first zone 
of 25, 50, 75, 100 and 200. As may be seen, the mobility of zone one has 
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FIGURE 14 Effect of Zone One Mobility Variation 
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little effect on the time to breakthrough, but does affect the amount of oil 
produced. For the streamline approaching the production well fastest, the 
conductivity ratio changes as the most advanced zone is produced. Equation 
4. 1 was expanded to include all zones considered in this paper. As may be 
seen, if the leading zone is produced, the conductivity ratio expression de- 
creases by one term in the denominator since it no longer exists in the 
reservoir for that particular streamline. Since the mobility of zone one is 
the same as the standard, the conductivity ratio approaches unity as each 
zone is produced, regardless of the value used as the mobility of zone one. 
Thus the mobility of zone 1 has little effect on the streamlines that travel 
the shortest distance to the production wells. The effect on the production 
may also be explained by Equation 4.1 in the same way as the stopping of 
the burn front was determined. However, the larger the mobility of the first 
zone, the larger the value of the denominator and the smaller the conductivity 
ratio. Thus, the greater the value of the mobility of zone 1 the shorter the 
distance from the injection well to the final position of the burn front. 

The mobility of the other zones may also change during the in-situ 
combustion process, but their effect should not be so significant as a change 
in zone 1 mobility since they are filled with a greater percentage of liquid 
and the temperature decreases rapidly between combustion front and zone 2. 

Injection Rate Variation 

Additional tests were conducted to determine the effect of increased 
injection rates on the total production and the time of breakthrough of the 
combustion front. For these tests, the relationship between injected air 
flux and the advancement of the combustion front was as developed by Martin. 
The results of these tests are shown on Figures 15 and 16. These results 



Total Production, bbls 
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Injected Air Flux, mcf/hr 




FIGURE 15 Effect on Oil Production by Variation in 
Injected Air Flux 



Time for Breakthrough of Combustion Front, Days 
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FIGURE 16 



Effect on Time of Breakthrough by Variation 
of Injected Air Flux 



41 



were rather interesting in that at higher injection rates the total production 
was lower. The explanation of this may be found in Figure 7. At high in- 
jected air flux rates the combustion front moves rapidly between the source 
and sink. The movement is much faster than the case of lower injection 
rates and results in an early breakthrough. As a result, the streamlines 
off the direct line between the injection and production wells are not able to 
move very far from the injection well. Thus, a smaller amount of oil is 
produced. The opposite action occurs at lower injection rates. The stream- 
lines on the direct line between the injection well and production well move 
rapidly, but do not breakthrough at the production well before the other 
streamlines have moved a considerable distance from the injection well. 

The graphic evidence of this action is shown in Figure 17 where the broken 
line represents the location of the combustion front for an injected air flux 
of 4000 scf/hr, and the solid line represents the combustion front for an 
injected air flux of 600 scf/hr at breakthrough. 

The results of the tests on the effect of increased injection rates in- 
dicate that economic considerations must be made to determine the maximum 
injection rate versus the time and amount of oil recovered. 

Bounded System 

Heretofore, the reservoir has been considered infinite in length and 
width. Now the reservoir will be examined with a boundary existing just 
outside the production wells. The pattern boundary is shown in Figure 18. 

In previous tests we noted the decrease in oil production due to the 
fingering effect and rapid ex^Dansion of the combustion front, thereby forcing 
oil out of the pattern area thus making it non-recoverable. As is shown on 
Figure 11, a majority of the oil within the pattern has been pushed away 
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— Injection Rate of 4, 000 ft^/hr 
Injection Rate of 600 ft^/hr 




FIGURE 17 Combustion Front Location at Breakthrough 
with Variation in Injected Air Flux 
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FIGURE 18 



Typical Pattern Showing Imaging Technique 
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from the injection well and in fact away from the production wells. In order 
to close the system, the imaging technique will be utilized. The imaging 
technique is a well-known method for confining streamlines inside a de- 
sired boundary. The technique is described in detail by Collins (13). In 
cases of constant reservoir thickness only the flow rate and location of each 
well is necessary. 

The pattern used, including the image wells is shown on Figure 18. 
The location of the burn front at breakthrough is shown on Figure 19. Com- 
paring the bounded system which has the same properties, injection and 
production rates as the basic model shows that the start of production and 
breakthrough occur fastest in the bounded system, and that the production 
rate of the bounded system lags that of the unbounded system. This could 
be expected since the imaging technique prevents flow across the established 
boundaries which are shown by the broken line in Figure 18. As a result, 
the oil laden zones are directed towards the production wells; however, in 
the unbounded system, the zones ahead of the burn take a more circuitous 
route as shown in Figure 1 and sweep part of the reservoir outside the area 
of the test. 

Reviewing the production data in Table in for the bounded system 
shows the start of production and breakthrough of the combustion front occur 
earlier than in the basic unbounded system. The bounded reservoir pro- 
duces about one-third as much oil as the unbounded. The importance of the 
existence of natural boundaries is therefore apparent. The existence of 
faults close to the pattern could cause similar reductions in production as 
the artificial boundary utilized in the above test. Fortunately, Lin (16) has 
developed a technique for construction of boundaries for the streamline 
model that may be used to simulate the shape of natural barriers found in 
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FIGURE 




o Injection Well 
o Production Well 



Location of Combustion Front in Bounded 
System 



TABLE m 



RECOVERY PERFORMANCE FOR BOUNDED SYSTEM 



Gas Produced 
(it^) 


Oil Produced 
(bbls) 


Water Produced 
(bbls) 


Time 

(hrs) 


0.0 


0.0 


0.0 


6.3 


0.0 


0.0 


0.0 


68.7 


19.0 


3.3 


3.3 


131.0 


113.0 


19.5 


19.7 


194.0 


223.0 


38.8 


39.5 


256.0 


334.0 


58.4 


59.7 


319.0 


453.0 


79.7 


81.6 


381.0 


556.0 


97.8 


100.0 


444.0 


668.0 


118.0 


121.0 


506.0 


767.0 


136.0 


139.0 


569.0 


858.0 


152.0 


156.0 


631.0 


930.0 


164.0 


169.0 


694.0 
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reservoirs. By including boundaries, the area swept and total production 
may be determined for various locations of the injection and production 
wells. 

Comparison with Field Tests 

To establish the ability of the streamline model to determine the 
location of the combustion front, data from two field trials were used as 
input data for the model. The first test utilized information provided by 
the Sun Oil Company (17). The test was conducted in the May-Libby reser- 
voir which has the characteristics as shown in Table IV . The pattern used 
was an inverted five spot with dimensions as shown on Figure 22. Of the 
information provided on Table IV, only the porosity, injection rate, and well 
pattern and spacing were utilized. The rationale for using only this infor- 
mation will be explained at the end of this section. The saturation and mo- 
bility of each zone was assumed to be the same as that of the basic test in 
the first section of this chapter. In addition the amount of reservoir hydro- 
carbon existing and consumed in the propagation of the combustion front was 
the same as in the basic test. The burned area is shown as the shaded 
area in Figure 22. As may be seen, the area burned in the field test and 
that indicated by the streamline model are extremely close. 

The second test was a comparison of the streamline model and a 
field test conducted by the Magnolia Petroleum Company (9). Table V 
shows the characteristics of the reservoir and Figure 23 shows the dimen- 
sions of the inverted five spot pattern. As in the previous test, the satur- 
ation and mobility of each zone were as used in the basic test as well as 
the amount of reservoir hydrocarbon existing and consumed. The burned 
zones are quite similar except in the region between the injection well and 



TABLE IV 



CHARACTERISTICS OF SUN OIL 
COMPANY TEST 



*Well Pattern: 


Inverted Five Spot 


*Well Spacing: 


As Shown on Figure 10 


* Average Porosity: 


31.2% 


Depth of Pay: 


3, 400 ft 


Thickness of Pay: 


4.4 ft 


Average Reservoir 
Permeability: 


1 , 069 md 


Average Amount of 
Oil in Place: 


740 stb/acre-ft 


* Average Injection 
Rate: 


8, 000 cu. ft/hr 



Used as input for streamline model 
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FIGURE 20 



Streamline Model Burned Zone versus 
Sun Oil Company Field Test (9) 
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TABLE V 



CHARACTERISTICS OF MAGNOLIA PETROLEUM 
COMPANY FIELD TEST 



*Well Pattern: 


Inverted Five Spot 


*Well Spacing: 


As Shown on Figure 23 


* Average Porosity: 


27.2% 


Depth of Pay: 


195 ft 


Thickness of Pay: 


17 ft 


Average Reservoir 
Permeability: 


7, 680 md 


Average Amount of 
Oil in Place: 


1, 350 bbls/acre-ft 


* Average Injection 
Rate : 


3, 200 scf/hr 



Used as input for streamline model 
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FIGURE 21 Streamline Model Burned Zone versus 
Magnolia Petroleum Company Field 
Test (17) 
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production wells land 4. The authors of the article published on the burn 
indicated that the reservoir properties dictated the non-radial formation of 
the burned zone. The streamline model of in-situ combustion is based on 
a uniform permeability; therefore, if the permeability of the areas in 
question are lower than the rest of the reservoir, little of the injected air 
would go in the direction of wells number 1 and 4. In addition, if the reser 
voir was tighter between the injection well and production wells 1 and 4 
the production of these two wells would be lower than that of the other wells. 
The authors confirmed this fact. 

A third test was attempted using the information published by the 
U.S. Bureau of Mines with regard to an in-situ combustion experiment in 
oil shale (11, 12). The result was not satisfactory because of the nonuniform 
ity of the permeability of the oil shale due to extensive fracturing around the 
injection and production wells. The test pattern consisted of an inverted 
five spot with the first set of four production wells approximately twenty- 
three feet away from the injection well. A second set of production wells 
approximately fifty feet away from the injection well were located on the 
same axis as the first set. With the extensive fracturing consisting of 
electrolinking, explosive fracturing and by hydraulic fracturing, areas of 
highly permeable oil shale existed and enabled excellent communications 
between the injection and production wells. Since the permeability in the 
rest of the reservoir was extremely low, all injected air would flow in the 
direction of the least resistance, i.e., toward the production wells. Making 
the assumption that the oil shale was a homogeneous, isotropic deposit, as 
was an assumption for development of the streamline model, was therefore 
erroneous. As a result, the test by the U.S. Bureau of Mines could not be 
simulated by the streamline model. 
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The articles relating the two field trials did not contain all the in- 
formation necessary to reproduce the field tests. As was stated above, the 
porosity, well spacing, and injection rate were the only data taken directly 
from the field tests. This was done for numerous reasons. First, the 
streamline model assumed that the reservoir was only one foot thick. Thus, 
the combustion front could be assumed to be a vertical front. Although this 
assumption is valid for a thin reservoir section, Moss has found that the 
assumption cannot be made for a thick reservoir (9). Secondly, the perme- 
ability is tied directly to the mobility of each zone. Since this area is one 
which is unknown at the present time, the permeability of the reservoir is 
of no use. Lastly, the difference in production, as demonstrated by the 
comparison of the Moss and Martin relationships between rate of advance 
of the combustion front and the injected air flux, would prevent meaningful 
comparison since it is not known if the aforementioned relationship is the 
same for all reservoir materials. 

Discussion 

The comparison of the streamline model to field trials demonstrated 
the model, as it exists, can approximate the location of the combustion 
front. Some of the variables that affect the combustion front location have 
been examined and the effect determined. The effect of other variables that 
warrant additional study are listed in the next chapter. 

As mentioned earlier, no attempt was made to determine the loca- 
tion of the oil bank or to duplicate production data due to the complexity of 
variables necessary to determine the velocities of zones 2 and 3. 



CHAPTER V 

CONCLUSIONS AND RECOMMENDATIONS 



Conclusions 

As a result of this study, the following conclusions may be 

drawn : 

1. The streamline model can be used to study the effects of vari- 
ables on the process of in-situ combustion. 

2. The streamline model can be used to determine the location of 
the burnfront with an accuracy dependent upon the accuracy of the input 
data and the homogenity of the reservoir. 

Recommendations 

1. The relationship between the rate of advance of the combustion 
front and the injected air flux should be determined for the specific reser- 
voir materials in order to increase the accuracy of the streamline model 
for in-situ combustion in the specific reservoir. 

2. A study of the problem of determining the saturation of various 
zones is necessary to obtain more realistic values of the zone composition 
and its mobility. 

3. Larger pattern areas should be utilized to determine the best 
economic relationship between injection rate and pattern area. 

4. An expanded study to include variation in reservoir thickness 
would be of great value . 
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5, The effect of dip should be considered as it affects the flow of 
oil during the combustion process, 

6. The results of this study were based on a reservoir with uniform 
permeability. With the extensive use of fracturing, a study comparing the 
effects of permeability changes is warranted. 
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Plot Location of Wells 
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APPENDIX B 



COMPUTER PROGRAM 



PROGfJAM COh'R ( !r.PUT»Ot)TPUT ) 

THIS HROGf!«M Pl.OTS iHt LOCATION Of- THE BURN FRONT 
dimension D (7 ) ,X.( ( 7) , Yw ( 7 ) tNSl. (41 f iTSUM (7 »4 1 ) iX(7*41) »Y(7i4i) , 
C0(7*A1) ,Ci-(7t4l) *C'W(7.41) ,CY(7 o4 1 ) .XEB(7«4l) »v£fi(7t4i) » 

TBi (7 »A1 ) ,XFd (4l ) »YFB (4 1) , W^^TPRU ( 7 ♦ 4 1 ) , XE2 ( 7 « 4 1 ) • YE2 < 7 . 

An «^F2 (41) I YF2 (‘♦I ) tCiASPPO (7,4 1 ) , XL 3 ( 7 , 4 1) , YE3 ( 7 , 4 1 ) , 

XF3(41) ,yF3( 41) ,AC4( 7,41) ,YL4 (7,41) ,XF 4(41) ,YF4( 41) 

READ 330, NST,NP,NI,NUAT,RI,FHP,pUR,ZM,PM 

READ 390, PVOLf'X,H,XHAX,YMAX, IPV 

READ 430, Sloped, BONE, fluxm 

READ 440, SLOPET,BTWO,i- LUXMN 

READ 5B0, NH0,FUEL,&ALP1 N 

read b?6, RATlo2»RATl03,RATl04,wRATl2,WRATl3,WRATl4,oRATI2,r,RATl3, 
1 gratu 

READ 5io, voL02,VOL03 

PRINT 500* NST,NP,NI,NUAT,Ri,Pcp,4v,PM 
PRINT 490, PV0LMX,H, XMAX,YKAX, IPV 
PRINT 450t SL0P£0,B0 n£ 

PRINT 530* SLOPtI,BTwO 
PRINT 4f>0, RH0,FUEL,GALPTN 
PRINT 470» RATt 02,RAT103,HATI04 
READ 540, SM8,TMB,F0RMb,FlVPB 
PRINT 350, FMB,SRB,TMb,FCRMR,FlVMb 
1T0T=NP*NI 

READ 340, <Q(I) ,XW(I) ,YW(I) ,1=1 ,iT0T) 

READ 540, Dmin,QM1NP,AREA,SdTST 
PRINT 570 

PRINT 3PCI* (Q(D ,XW(I) ,YW(I) ,lsi,lT0T) 

PRINT 460 , QMlN,(iMlNP, AREA,SDTST 
SVALPV = AREA«HopOK/QMIN,»3600, 

SVALUE = IPV‘'SVALPV/fJOAT 

correction factors from English to metric units 

6ASPPD(j,K)=0,0 

WAIPPD ( j ,K) =0.0 

M1)STX = MDSTX«30.48 

HDSTY = MCSTY,*30.4cI 

XMAX = XmaX<^30 .48 

YMAX=YmAX''^30.48 

SDTST = SCTS ( *>30, 48 

H=H«30,48 

RI=SCTST/50.0 

NTOTsNI ♦^Jp 
PO=50.«-30.48 
OMlNiOMir)*> ^,88o55 
0(1INP = OMINPo 7. 88055 
NIl=NP*l 

CTR=1.0 

PI=0.0 
f'P = 0,0 

DO lO 1=1,IT0T 

correction factors fpom English to metric units 

XW ( I ) =X\v ( n *>30 , 48 
YW(I)=YV,(in>30.48 
0(1) ={v'(I)0/. 88055 
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f'I=Pl” <0 ( P <^ZM) / ( <''AL0G ( ( Xw ( N 1 1 ) *R 1 - X ». ( I > ) **?♦ 

1 ( YW (NI 1 ) -YW ( 1 ) ) 

PP = PP- (0 ( 1 ) *»ZM) / (12,5664«FMB0H0PM) «AL0G ((Xw(l)*RI-XW(I))<*«2* 

1 (YWU)-Yw(I) ) 

10 CONTINUE 

CALL WGNPLT 

CALL PLT ( i0.0t5,5*-3> 

CALL PLT (-3.0,~3.0t3) 

CALL PLT (-3.0,3.0i2) 

CALL PLT (3, 0.3. 0.2) 

CALL PLT (3. 0.-3. 0.2) 

CALL PLT (“3. 0,-3. 0.2) 

SCL=, 0010024783 
DO 20 1 = 1. noT 

CALL GYPBOL ( Xw ( I ) «SCL « YW ( 1 ) USCL . i 2 1 . 1 . 0 . . - 1 ) 

20 CONTINUE 

DO 40 JrN'Il.NToT 

NSL (J) =NST*Q( J) /UM IN. 0.45 

PRINT 370 

1B=NSL ( J) 

DO 30 K=1.IB 

X(J,k)=RI*COS(1*6.2832<^K/NSL( J) 

Y (U.K) =Rp'5lN( l.b.2832«K/NSL ( J) ) *YW ( J) 

XE2( J.K)=X(J,K) 

YE2 ( J.K) = Y U.K) 

XE3 ( J.K) =X ( J,K) 

YE3 ( J.K) =Y (J.K) 

XE4 ( J.K) =X ( J,K) 

YE4 ( J.K) =Y < J.K) 

OW( J,K)r0.0 
CO(J,K)=0.0 
CW(J,K)=o.O 
TBT (J.K) =0.0 
TSUM( J.K) =0.0 
30 CONTINUE 
40 CONTINUE 
50 CONTINUE 

DO 200 J=NI1.NtOT 
lB = r;<:L ( j) 

DO 280 K=1.IQ 

IF (GW(JtK) .GT.O.O) GO TO 260 
XFBrX ( J,K> 

YFB=Y ( J,K> 

XF2=XE2( J.K) 

YF2=yE2( J.K) 

XF3=XE3(J.K) 

Yf3=YE3( J.K) 

XF4=xL4(J,K) 

YF4=yL4 ( J.K) 

CALL PLT (X(J,k)"SCL,Y( J,K)*SCL.3> 

TSUMI=TSUM( J.K) 

IF (AbS (XFB) .GT.XmAX) 270.60 
60 CONTINUE 

IF (ABS (YFB) .GT.YmAX) 270.70 
70 continue 
VELXP=0.0 
VELYP=0.0 
VELX?=0.0 
VELY2--0.0 
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VELX3=o.O 

VELY3=0.0 

VELX4=0,0 

VELY4=0.6 

PB=0,0 

P2=0,0 

P3=0,0 

P4=0,0 

PK5T = /:m/ ( lii,5664«FMn*H"PB) 

VK5T = ZM/'(6t2B32*Ht’P0R<*PN') 

(jO 80 L = l»iTOT 

FAC4=(XK4-'^W(L) )*»2* (Yf-4-YW(L) 

P4 = P4-PkST*Q (L) <*4106 (F AC^> 

VELY4 = V(-LY4*((vKST'»Q(L)/CTP)O(YF4-YW(U)/FAC4)/P0P 
VELX4-:velX^* t < VI^ST^Q (L) /CTR) ** (XF4”XW (L) ) /FaC4 ) /POR 
VELS4 = SCP1 < V'ELX4<‘'^2<-VE:LYA<»«2) 

FAC3=(XF3“XW (L) (YF3“Yw (L) )««2 

P3 = P3"PkST<‘u (L) <>ALOG (FAC3) 

VELX3=VFLXJ+ ( ( VXSTOQ (L) /CTR) * (XF3“Xw (L) )/FAC3) 

VELY3 = Vri < (VXST*Q(U/CTR)o(YF3"YW(L))/FAC3) 

VELS3 = 5GRT < VELX3<‘o2*vELY3t "2) 

FaC 2 = (XF2-XW (L) ) ( Yt- 2-Y/; (L) ) 

P2 = P2“PKST<*U (L) <-'ALOG (F ac2 ) 

VELX2 = VELX2* ( ( vxbT <»Q ( L ) /CT R ) <» ( X F?-XV/ (L) )/FAC2) 

VELY2 = VEL Y2* ( ( vKbT^Q (L ) /CTP) " ( YF2“YW (L) ) /FAC2) 

VELS?--srRT ( VELX2««2 *vLLY2“*2) 

F ACB-- ( XF0“XVML) ) ( Yr 8-Yw (L) ) t»o2 

PR = PP-PKST‘‘0 (L ) *ALOG (F ACB ) 

VELXP = VELXb* < (vKST«Q (L) /CTR) * (XFR"XW (L> ) /FACB) 

VELYP=VELYB* < (VKST**0 (L) /CTR) *> < YFB-YVv (L) ) /F'ACB) 

60 COMINUE 

VELSP = SGRT (VELXB<»«2*vELyB<J02) 

IF (SQR] ( (AHS(XF4) -PQ)»*2* (APS(YF4)-PQ)«<»2) .LT.Pl ) GO TO qO 
CTY= (Pp-pl ) / (Fpb" ( ( (PB-PI )/FMR) ♦ ( (P2 “PB) /SpB) ♦(<P 3-P2)/TyB)*(<P4- 
1 P3) /F OP.MB) ♦ ( (PP-P4 ) /F I VMp ) ) ) 

60 TO 120 
90 COMINUE 

IF <SOP!( (AHS(XF3)"PQ)<‘»2*(ARS(YF3)-PQ)«*2) .LT.RI) GO TO lOn 
CTY = (PP-pl ) / (Fpfi« ( ( (pb"Pl ) /FM8) ♦((F2-PH)/SyB) + ((P3"P2)/TPB)*((PP- 
1 P3 )/FoRM8))) 

GO TO 120 
loo COKTlRUE 

IF (SORT ( (AbS(XF 2)-P0)««2* (APS(YF2) -PC)o<»2) .LT.RI) GO TO 1 1 6 
CTY= (Pp-pl ) / (FmB*» ( ( (PB“P1 ) /FM8) ♦ ( (P2-PB) /SpB) ♦ ( (P 3-'P2 > / TPB ) ) ) 

60 TO 120 
110 CONTINUE 

UEM = FMR«» ( ( (PP-pB) /FM8) ♦ ( (P0-PI) /SF’S) ) 

CTY= (Pp-pl) /OEu 
120 CONTINUE 

VELSA=VElSB«PM/2M«CTY 

IF (VELbA.LT»FLUXM) 140il30 
130 CONTINUE 

velb=slope.o«velsa*bone 

60 TO 1530 
140 CONTINUE 

IF (SLOPETovELSA.lE.ABS(BTwC) ) GO TO 270 
VELB=SLOPtTovELSA*BTwO 
IbO CONTINUE 

AMTO 1L = G ALP TN«RH0«, 00006684 
POR=POR-aMU)IL 




M 
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FACT 2 = Af^T0lL‘‘(l. 0-FUEL) « < V0L02/ ( VOL 02 + V 0 L 03 ) ) / (PnR«vnL02) 

FACT3 = APT01L* ( 1 ,0-FUEL) (V0L03/ ( v0l 02 ♦ V0L03 ) ) / (P0P"vnL03> *FaCT? 
F>0H=,i:5 

VELFACrVELW/VELSa 
SVELXArVEL^A 
SVELY‘» = VELY4 
VELTA=VELS'» 

IF (?ORT( (ARS(XF3)-PQ)«*2* (ARS<YF3)-PQ)«"2) .LT.RI) GO TO 16o 
SVELX3 =v£L^3*''VELF AC<>FACT3/poh 
SVELY3=V£LT3'-»V£I.FAC<-'FACT3/P0R 
VELT3 = SCRT <SVEL)^3«*'2*SVELY3‘>*‘2) 

160 CONTINUE 

IF (SORf ( (AbS(XF 2)-PO)<‘<»2* (APS‘YF2)-PO)‘»«2) .LT.PI) GO TO l7o 
SVELX2 = vELa2'WELFAC«’FACT3/POR 
SVELy 2=VELY2<»VELF AC«FAC1 2/POR 
VELT2 = S0RT (SVELX2<“>2*bvELY2<**‘2) 

170 continue 

SVELxB=vEL''^BoVELF AC/POR 
SVELYR = VELYH<>VELFAC/pOR 
VELT()=SGRT (SVEL^H<**2^i>VELYB«’<»2) 

VELTr = VELT<t 

IF (vElTA.OT.VeLTM) VcLTNrVELTA 
IF (VELT3.^T, VeLTM) v£LTI<' = VELT3 
IF (VELT2,OT,VeUM) vELTB = VELT2 
IF (VElTB,LT.VfL1m) vELTF=VELTP 
IF (PB.^T.O.y) GO TO 180 
DO 180 L=1»UP 

RAD = SURT ( (XFB-X.V (L) ) <*^2* (YFB-Yw (L> )®«2) 

IF (RAO.GT.RI) go to 180 
TBT (J,K)=TGUMI 
OW( JiK)=qmIn/NSL( J) 

GO TO 2?0 
180 CONTINUE 

ut=ri/veltm 

IF (TSU^^♦OT,GT.SVALUE) 190. 20<) 

190 CONTINUE 

DT=SvALUE“iSUMI 

200 continue 

XF0=xFR,DT^SVELXa 
YFB=yFR*dT*SVElYB 
XF2=XF2*nT"SVELX2 
YF2 =yF2*0T^SVElY2 
XF3=xF3*OT"SVElX3 
YF3 = YF3*DT**SVElY3 
XFA=xF4*dT*SV£lX4 
YF4 = yFa*0T*'SVElYA 
CALL PlT (XFB^^SCL.YFB^SCL.E) 

TSUMI=TSUUI *DT 
SOTST=70»7"30,40 

IF (S’JRH (ABS(XF4)-PQ)««2+ (ABS<YFA)-PQ)«<»2) .LT.RI) GO TO ?l0 
RATIC=0,0 
WRAT I0=o. O 
GRATlO=o ,0 
GO TO 240 
210 CONTINUE 

WRATlOrWRAl I4 
RAT10=RATI04 
GRAT i0=GRAT 14 

IF (SORf ( (A rS(XF3)-PQ)««2* (ABS<YF3)-P0)««2) .LT.RI) GO TO 220 
GO TO 2‘*0 
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2H0 C0HTI‘'^UE 

RATlC=PflTl^3 
wRATlOrwRAT 13 
ORATlOrGRAT 13 

IF (SQRI ( (ABS(XF2)“PQ) (ABS ‘YF2)-PG) »*2) .LT ,RI ) GO TO ?3o 
GO TO 2^0 
230 continue 

RAT IO'-RAT 1^2 
WRATlOsWRAl 12 
GRAT lOnGRAT 12 
2A0 CONTINUE 

AODPRO = OT* UjMIu* (ZM/pM)/ (NSL ( J) oVtLSB^POR) I^I.OOOOOlRl^lVELSR'^PORo 
1 CTY*^ tPM/Zf-'') ) *.0003^9) 

CO(J,K)=COU»K) +ADDPR0‘*RATI0 
GASPRDI J,K) =GASPR0( J,K) ♦ ADOPROoGpAT I 0 
KATPR0( J,K) = WATPr(0 ( JtK) tAODPROowRAT 10 
IF (ABS(XFB) .GT.XmAX) 270*2=0 
250 continue 

IF (ABS (YFB) ,GT. YMAX) 270*260 
260 CONTINUE 

IF (TSUP'I .Lt.SVALUE) go TO 70 
CALL SYPBOL (XFd<^SCL*YFH«SCL*. 07*11*0, *-1) 

270 CONTINUE 

X ( J*K) =XFB 
Y(J»k)=YfB 
XE2( J*K) =XF2 
YE2( J*K)rYl-2 
XE3(j*K)=xF3 
YE3( j*K)=Yf-3 
XEA(J«K)=XfA 
YEA(J*K)=Yi-4 
TSUM( J,K) =TSUKI 
CY(J*K)=CTY 

2B0 continue 
290 CONTINUE 
COP=0«0 
C6P=0,0 
CWTRP=0,6 
CTR=0»0 
CwP=o»0 

00 3C0 J=Nll*NTOT 
IB=NSL(J) 

UO 300 K = 1 * IB 

CWP = CWP,QW ( J*K) (SVALUE-TBT ( J*K) )*. 0000062898 
C0P=C0P*C0 (J*K) <>.0000062898 
CGP = CGP,GABPRD(J*K) ‘►.000035314 
CwTRP=CWTRP*WAtPHD (J*N) 0000062898 
CTR=CTR*CY ( J*K)/1B ' • 

300 CONTINUE 

TINE=SVALUt/36oO.O 

CLP=COP,CwP 

PRINT 360* CGP,C0P*CwTRP*T 1ME 
IF (CWPiGT.0,0) 320*310 
310 CONTINUE 

SVALU£ = SVALUE*SVALPV/Nt)AT*10,0 
60 TO 50 
320 CONTINUE 
PRINT 550 



.'."T " - 
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PRINT AOO* NST,NP*NI ,NUAl tRi ZM,PM 

PRINT 5f0 

PRINT 410* PVOLMX, H, XMAX* YMAXf IPV 
PRINT 570 

print 420* (Q(D »XW(1) *YI«m *I = l,ITOT) 

CALL End PL T 



330 FCHMaT (4i5*4F5.2* 

340 FORMAT <F'io.O*5X*F 

3b 0 format ( 16X,«;-0H1LI 1 Y 

1 16Xf«>MU!3ILlT Y 

2 l6Xt<‘H0BILlTY 

3 l6X,<vMO(3ILIlY 

4 16X**M03ILITY 



F 7.0) 

5,0«BX,fB«0) 

CF FLUIP IK ZONE 1 

CF FLUID IK ZONE 2 

CF FLUIU iK ZONE 3 

CF FLUID Ik ^0N£ 4 

CF FLUID Ik ZONE 5 



‘^FlO.B*/* 
*F10.3*Z* 
*F10.3./« 
“Flo. 3*/* 
«F10.3) 



360 


F ORMAT 


(/16X,El0.2,bX,E10.2,ex.Ll0.2*12X,£l0.2) 






370 


format 


(///16X.OGAS RRCDUCED OIL PRODUCED 


WATER PRoDUC 




lEO 


time IN H0URS«*Z* 16X,**cOBIC FEET 


barrels 






2 Barrels^) 






3b0 


format 


( /16X,F10.1.5X,F10il»5X.F5.1) 






390 


format 


(4F7.2* 17) 






4U0 


format 


( /i6X,4l5*10E10,3) 






410 


format 


(/16X,5E10.3) 






4ii0 


format 


( /16X,e10.3*5X,£10,3,5X,E10.3) 






430 


format 


(3F10.8) 






440 


format 


(3F10.8) 






4b0 


format 


(//16X»17HSl0PE like Oke = F9 . 7 ♦ 3X i 1 4 HL InE CONST, = 


FI 




lu.8) 








460 


format 


(16X,«MINIMUM INJECTION RATE = *El0.2*/* 


16X»<>MINIMIJM 


PRO 



IDUCTiOn rate = *EI0.2*/» 16Xf«wELL FATTERK AREA = «F 1 0 . 2 . / * 1 4 X ♦ 

• 2 ^distance between INJECTION WELL AND PRODUCTION WFLL = “FI0.2) 

470 format (16X»«^AM0UNT OF OIL IN ZONE 2 = «F 5 . 3 • / * 1 6 x ♦ AMOUNT cF 0 
IIL IK zone 3 = «F5.3»/* IfcX.oAMCuNT OF OIL IN ZONE A = «F5.3) 

4b0 format (/16x*«AVERAGE density UF RES.(RhO) = oF10.4»/,16x.<*A'-' 0U 
INT OF OIL CONSUMED AS T UEL = <>F 1 0l A , / , 1 6X * «> AVER AGE RpS. OIL COnTEN 
2T GAL. per' TON = «F10.4) 

490 FORMAT * ( Zl 6X * <>RESER VO I R PORE VOLUME = «F9 . 2 * * / * 1 6 x * »TR ICkNE «:S 0 
IF reservoir = oF6.2,/*16X,*^MAXIMu'^ LIST IN X DIRECTION = <>Fa.2,Z. 

2 16X.*JMAXIMUM UlST ]N Y DIRECTtO'^ = *F6 . 2 * / ♦ 1 6X i « I K jE CTED PORP VO 
3LUMES = «ID 

500 format (/I6X,i‘NUMBER OF STREAM TUBES = » 1 5 * / ♦ 1 6X * «K'UyBER OF PROD 
lUCTlCN KELLS = *■16*/. IGX.ONUMHER OF INJECTION wElLS = I 5 * / . 1 6X * 
2DAT :: «>I5*Z*16X.<:RADlUb CF WELL BORE IN FEET = <>F5 . 2 . / * 1 6X . opQroS I 
3TY OF RESERVOIR = <»r5. 2 . / * 16X , ocuMPRESS I B IL IT Y FACTOR = “F5.2./.16 



4X, ^reservoir pressure in atmospheres = <»F7.0) 



510 

520 

530 


F ORMAT 
FORMAT 
F ORMaT 


(2F10.4) 

(9F5.3) 

(//16X*17rtbLCPE LIKE = 


F9.7*3X.14HLImE 


CONST. 


= FI 


540 

5b0 


K.8) 

F ORMAT 

format 


(4F10.3) 

(/18X*«NST KP NI NCAT 


RI Fm0 


POP 


RAT 1 


560 


10 ZM 

format 


PM«) 

(//1ox*«pvoLmx h xmax 


YMAX IPV 


AM8«) 




570 


format 


(//16X*» well rate 


X LOC Y lhC 


«) 




5B0 


F OHMAT 

end 


(3F10.4) 









APPENDIX C 



NOMENCLATURE 

2 

A Cross sectional area, L 

C Constant 

F Factor 

h Thickness, L 

2 

k Permeability, L 

H Length, L 

n Moles, dimensionless 

9 

P Pressure, M/LT 

3 

Q Injected volume, L 

q Injection rate, L'^/T 

2 2 

R Universal gas constant (per mole), ML /T t 

SL Slope constant 

S Saturation, dimensionless 

T Temperature, t 

t Time, T 

u Flux, L/T 

V Volume, L^ 

V Velocity, L/T 

z Compressibility factor, LT^/M 

Greek Letters 

y Conductivity ratio, dimensionless 

X Fluid mobility, L^T/M 
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II Viscosity, M/LT 

$ Potential, M/L^T^ 

(j> Porosity, dimensionless 

77 Constant = 3. 1416 

p Density, M/L^ 

Subscripts 
B Bulk 

b Burn 

eff Effective 

g Gas 

1 Counter for number of zones 

k Counter for number of zones 

m Average 

n Maximum number of zones 

res Value at reservoir conditions 

s Standard fluid 

sc Value at standard conditions 

T Total 

w Water 

X x-coordinate 

y y- coordinate 

z z- coordinate 

2 Zone 2 

3 Zone 3 

4 Zone 4 

5 Zone 5 
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ADDRPO 


Total fluid accimiulation for streamline under inves- 
tigation, L'^/T 


AMTOIL 


O 

Amount of oil in streamline under investigation, 


AREA 


Surface area bounded by production wells, 


BONE 


Value of C in Equation 3.8 


BTWO 


Value of C between FLUXM and FLUXMN 


CO 


Cumulative oil produced by a streamline, 


COP 


Cumulative oil produced, field-wide, 


CTR 


Pattern conductivity ratio for producing well, 
dimensionless 


CTY 


Conductivity of individual flow tube, dimensionless 


CW 


Cumulative production of injected fluid 


CWP 


Cumulative production of injected fluid field-wide 
basis, l3 


CY 


Counter for the conductivity of individual flow tubes 


DT 


Time required to move one step on a streamline, T 


FACB, 2, 3, 4 


Denominator within the summation of the velocity 
equations for each zone, l2 



FACT2, FACTS Term computed by Equations 3.37 and 3.41, 

dimensionless 



FIVMB 


Mobility of fluid in zone 5, L^T/M 


FLUXM 


Transition point of relation between air flux and 
advancement rate of burn, L^/T 


FLUXMN 


Minimum air flux required to move burn, L^/T 


FMB 


Mobility of injected fluid, L^T/M 


FORMB 


Mobility of fluid in zone 4, L^T/M 


FUEL 


Amount of oil in place that is consumed by the burn 


GALPTN 


Average gallons of oil per ton of reservoir material, 
L^/M 
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GASPRD 


Cumulative gas produced on a field-wide basis, L^/T 


GRATIO 


Amount of gas in zone at production well 


GRATI2, 3, 4 


Amount of gas in each zone 


H 


Thickness of reservoir, L 


I 


Counter for number of wells 


IB 


Counter for desired streamlines 


IPV 


Number of pore volumes injected 


ITOT 


Number of image wells 


J 


Step counter along streamline 


K 


Streamline counter 


L 


Counter for number of wells 


MULT 

NDAT 


Number of even increments of SVALPV at which 
production history will be accumulated 


NI 


Number of injection wells 


Nil 


First injection well 


NP 


Number of production wells 


NSL 


Number of streamlines emanating from a well 


NST 


Number of streamtubes emanating from a well 


NTOT 


Total number of real wells 


PB, 2, 3, 4 


2 

Pressure of leading edge of each zone, M/LT 


PI 


Pressure on well bore of injection well, M/LT'^ 


PKST 


Constant preceding summation in Equation 3. 5 


PM 


2 

Average reservoir pressure, M/LT 


POR 


Reservoir porosity, dimensionless 


PP 


Pressure on well bore of a production well, M/LT'^ 


PVOLMX 


Reservoir pore volume 


Q 


Volume of fluid, 
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QMIN 

QMINP 

QW 

RAD 

RATIO 

RATI02, 3, 4 

RHO 

RI 

SCL 

SDTST 

SLOPED 

SLOPET 

SMB 

SVALPV 

SVALUE 

SVELXB, 2 , ' 

SVELYB, 2, ; 

TBT 

TIME 

TMB 

TSUM 

TSUMI 

VELB 

VELFAC 



o 

Minimum injection rale, L'^/T 

O 

Minimum production rate, L /T 

3 

Volume of injected fluid, L 

Distance between point of interest and production 
well, L 

Portion of fluid in zone at production well 
Portion of fluid in zones 2, 3, 4 

O 

Density of reservoir, M/L 

Radius of well bore, L 

Scaling factor to reduce scale of well 

Distance between production well and injection well 

Value of SL in Equation 3.8 

Value of SL when air flux is between FLUXM and 
FLUXMN 

Mobility of fluid in zone 2 
(Pattern area)*(H)*(POR) /QMIN, T 
(IPV)*(SVALPV)/NDAT, T 

, 4 Velocity of each zone in x-coordinate direction, L/T 
, 4 Velocity of each zone in y-coordinate direction, L/T 
Time of breakthrough of driving fluid, T 
Time since start of injection of recorded production, T 
Mobility of fluid in zone 3, L^T/M 
Cumulative time since start of injection, T 
Counter for cumulative time since injection, T 
Velocity of burn, L/T 

Ratio of velocity of burn to velocity of air across 
the burn 

Air flux, L/T 



VELSA 
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VELSAR 
VELSB, 2, 3, 4 

VELTB, 2, 3, 4 

VELTMX 
VELXB, 2, 3, 4 
VELYB, 2, 3, 4 
VKST 

V0L02, 3 
WATPRD 
WRATIO 
WRATI2, 3, 4 
X 

XEB, 2, 3, 4 

XFB, 2, 3, 4 

XMAX 

XW 

Y 

YEB, 2, 3, 4 

YFB, 2, 3, 4 

YMAX 

YW 

ZM 



Velocity of the air across the front of each zone along 
the streamline, L/T 

Velocity of the leading edge of each front along the 
streamline, L/T 

The maximum velocity of the above four, l/T 

Velocity of each zone along the x-coordinate, L/T 

Velocity of each zone along the y-coordinate, L/T 

Constant preceding summation in Equation 3.6 and 
3.7, LVT 

Amount of oil in zone 2 and 3 

3 

Cumulative water produced, L 

Amount of water in zone at production well 

Portion of zone 2, 3, 4 that is water 

X-coordinate of fluid interface along streamline, L 

X-coordinate of each zone along streamline L that 
is to be stored, L 

X-coordinate of each zone along streamline L to be 
used in computations, L 

X-coordinate of maximum pomt of interest, L 

X-coordinate of well locations, L 

Y-coordinate of streamline L on well bore, L 

Y-coordinate of each zone along streamline L to 
be stored, L 

Y-coordinate of each zone along streamline L to 
be used in computation, L 

Y-coordinate of maximum point of interest, L 
Y-coordinate of well locations, L 

2 

Compressibility factor of driving fluid, LT /M 



BIBLIOGRAPHY 



1, American Petroleum Institute, Petroleurn Facts and Fijmres , 197 1 

Edition, Port City Press, Licorporated, Baltimore, MarylancI^ 

May 1971, p, 115, 

2, National Petroleum Council, U , S, Energy Outlook. National Petro- 

leum Council 1971, pp. 155-182, 

3, LeBlanc, J,L, and B, H, Caudle, "A Streamline Model for Second- 

ary Recovery”, Soc, Pet, Eng, J, Vol, 11, 7 (1971), 

4, Collins, R,E,, Flow of Fluids Through Porous Materials, Re inhold 

Publishing Corp, , New York, 19 d1. 

5, Rust, C,B,, A Streamline Simulator for Oil Bank Build-up in a 

Waterflood, M.S. Thesis, The tihiversit'y of Texas at Auslin, 1972, 

6, Howard F,A, , "Method of Operating Oil Well” U,S, Patent 

1, 47o, 348 to Standard Development Company, November 1923, 

7, Emery, L,W, , "Results from a Multi-Well Thermal -Recovery Test 

in Southeastern Kansas”, J, Pet, Tech, Vol, _12, (June 1962), 
pp, 671-678, 

8, Sheinman, A,B,, K,K, Dubrovai, N,A, Sorokin, M,M, Charvigin, 

S, L, Zaks and K, F, Zinchenlco, "Gasification of Crude Oil in 
Reservoir Sands”, Pet, Engr,, Dec, 10-27, (1938), 

9, Moss, J,T,, P,D, White and J,S, McNiel, Jr,, 'In -Situ Combustion 

Process — Results of a Five Well Field Experiment in Southern 
Oklahoma”, Trans, AIME (1959) 2 ]^, 55, 

10, Gates, C,F, and H, J, Ramey, Jr,, "Field Results of South Belridge 

Thermal Recovery Experiment", Trans, AIME (1958) 213, 235-244, 

11, Carpenter, H, C,, E,L, Burwell andH,W, Sohns, "Evaluation of an 

In-Situ Retorting Experiment in Green River Oil Shale", J, Pet, 
Tech, (January 1972), 21, 

12, Burwell, E,L,, T,E, Sterner, andH,C, Carpenter, "Shale Oil Re- 

covery by In-Situ Retorting — A Pilot Study", J, Pet, Tech, 22, 
(December 1970), 1521, 

13, Allred, V,D, , "Some Characteristic Properties of Colorado Oil 

Shale Wliich May Influence In-Situ Processing”, Quarterly of the 
Colorado School of Mines, ^ (July 1964), 



74 



75 



14. Martin, W. L., J.D. Alexander and J. N. Dew, "Process Variables 

of In -Situ Combustion", Trans. AIME (1958) 213, 28. 

15. Caudle, B.H., Fundamentals of Reservoir Engineerins , Soc. Pet. 

Engrs., AIME, Dallas, Texas, (1967). 

16. Lin, Jer-Kuan, An Image Well Method for Bounding Arbitrary Reser- 

voir Shapes ii TfTie Streamline Model. Ph.D. Dissertation, The 
University of Texas at Austin, 1972. 

17. Hardy, W.C., P.B. Fletcher, J.C. Shepard, E.W. Dittman and 

D. W. Zadow, "In-Situ Combustion in a Thin Reservoir Contain- 
ing High-Gravity Oil", J. Pet. Tech. (February 1972), 199. 

18. Merrick, R.J., Streamline Flow Solutions for Predicting Recoveries 

by Cycling Multiwell . Anisotropic. Stratified Gas Fields . Ph.D. 
Dissertation, The University of Texas at Austin, 1969. 



VITA 



Talbot Wade Bone was born at Montgomery, Alabama, on March 
21, 1944, the son of Charles W. and Dorthy J. Bone. He attended public 
school in Hummelston, Emmas, and Harrisburg, Pennsylvania, graduat- 
ing from Central Dauphin High School in 1962. He received his Bachelor 
of Science in Civil Engineering from North Carolina State University 
in 1966. 

On the 18th of March 1963 he enlisted in the U.S. Navy Reserve 
and in June 1966 received his commission as an Ensign in the United 
States Navy Civil Engineer Corps. After attending the Civil Engineer 
Corps Officer School at Port Hueneme, California, he was assigned to 
the U.S. Navy Mobile Construction Battalion Seventy-One as Assistant 
Company Commander and as a Company Commander. Subsequently he 
served tours as the Army/Air Force Construction Program Coordinator 
for the Officer In Charge of Construction, Thailand, and most recently 
as the Assistant Public Works Officer for the Naval Amphibious Base 
Little Creek, Norfolk, Virginia. 

He was selected for post graduate education in Petroleum Engi- 
neering and entered the Graduate School of the University of Texas in 
June 1972. 

Bone married Christina Van Horn of Harrisburg, Pennsylvania 
on September 6, 1964. They have two children, Todd, age 7, and 
Clark, age 4. 

Permanent Address; 3404 Derry Street 

Harrisburg, Pennsylvania 

This thesis was typed by Judith A. Schluter. 






I 



4 






I 



143009 

Thes I s 
B6794 Bone 

Streamline model of 
in-situ combustion in 
hydrocarbon bearing 
jureservoi rs . 



” ^ 143009 

Thesis 
B6794 Bone 

Streamline model of 
in-situ combustion in 
hydrocarbon bearing 
reservoirs. 



thesB6794 

Streamline mode! of in-situ combustion i 




3 2768 001 01718 9 
DUDLEY KNOX LIBRARY 









